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We present preliminary results about topology and several other quantities at non-zero baryon density and 
finite temperature in full QCD with N c = 2. 



1. INTRODUCTION 

The behaviour of the QCD vacuum at finite 
density can be tested in the heavy ion colliders. 
Therefore it is interesting to have a set of predic- 
tions from QCD and the lattice is a good means 
to obtain them in the non-perturbative regime. 

We have analysed the model with two colours 
in order to make use of standard simulation al- 
gorithms. We expect that this simplification may 
leave at least some of the essential features of the 
problem. 

We give a progress report of a simulation per- 
formed on a 14 3 x 6 lattice by using 8 flavours 
of staggered fermions at masses am = 0.05 and 
am = 0.07 at several values of (3 = A/gfi where 
go is the lattice bare coupling. The Monte Carlo 
algorithm is the standard hybrid molecular dy- 
namics [flf . The action is [01 
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\ V4(n){e a ^(n)U 4 (n)^(n+Z) 

n 

- C - a 'V(n + 4)[/ 4 (ri + 4)V(n)} 

+spatial terms, (1) 

where the "spatial terms" are those independent 
of the chemical potential fj, and 774 (n) are the stag- 
gered phases. 

Our results are preliminary because i) we have 
not yet calculated all the renormalization con- 
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stants necessary for the evaluation of the physical 
topological susceptibility x H an d ii) we have not 
calculated the physical units by determining the 
lattice spacing a. 

2. PARTICLE DENSITY 

In Fig. 1 we show the particle density p as a 
function of the chemical potential /1 for masses 
am — 0.05 and am — 0.07. In both cases we see 
a cubic dependence for large values of fj, that is in 
agreement with the analytical expectation [|J, 
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for /i T where T is the temperature. 

3. TOPOLOGICAL SUSCEPTIBILITY 

The topological susceptibility x is an important 
parameter of the QCD vacuum. In the continuum 
and for the SU (2) gauge group it is defined as 



\ = / d 4 x d M (0\T{K M (x)Q(0)}\0) , 
where K^x) is the Chern current 



KM = 



g is the QCD coupling constant and Q(x 
d^K^{x) is the topological charge density, 
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Figure 1. Baryon density p as a function of the 
chemical potential at (3 = 1.5. 



Eq. (||) uniquely defines the prescription for the 
singularity of the time ordered product when 

On the lattice, the evaluation of \ does not fol- 
low the above prescription. As a consequence we 
have to subtract the singularity which shows up 
as an additive renormalization. This subtraction 
is performed following the method of Rcf . (|] . 

It is well-known that \ undergoes a sharp drop 
at the deconfinement temperature, both for three 
colours H and for two colours || . We have mea- 
sured this quantity as a function of the chemical 
potential in order to study its behaviour at finite 
density as well as finite temperature. In Fig. 2 
we show that the behaviour is qualitatively anal- 
ogous to the above-described one at finite tem- 
perature and zero density. The quantity Z in the 
vertical axis is the finite renormalization which re- 
lates the topological charge in the continuum and 
on the lattice . We are presently running more 
statistics to sweep the sector beyond the critical 
density and to improve the error bars. 

4. THE POLYAKOV LOOP 

To investigate the deconfinement properties 
of the theory across the phase transition at fi- 
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Figure 2. Topological susceptibility as a function 
of the chemical potential at a finite temperature, 
(3 = 1.5 and am — 0.07. 



nite chemical potential, we have measured the 
Polyakov loop II. Our results, shown in Fig. 3, 
indicate that the theory becomes deconfining be- 
yond a critical density and suggest that the tran- 
sition occurs at afi ~ 0.4 for fixed temperature. 
A similar trend can be observed by varying the 
temperature at zero chemical potential, Fig. 4, 
where we have traded the temperature T for the 
lattice coupling (3. 



5. CHIRAL TRANSITION 

The chiral symmetry is known to be restored 
at the critical temperature. This result calls for 
a study of the behaviour of this symmetry across 
the critical density. In Fig. 5 we show the chi- 
ral condensate (ipip) as a function of the chemical 
potential. We see that there is a transition ap- 
proximately at the same value of \x at which the 
deconfinement took place. An analogous figure 
for (ipip) is obtained at zero chemical potential as 
a function of T. 
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Figure 3. Polyakov loop II as a function of the 
chemical potential for j3 = 1.5. 



6. PARTICLE SUSCEPTIBILITIES 

The particle susceptibility Xparticic is defined as 

Xparticic = ^2 G W ' ( 6 ) 
t 

where G(t) is the propagator of such a particle 
at time t. If for example we consider a pion 
and a single pole dominates the propagator, then 
~ Z^/rn^ where Z„ is the pion field renor- 
malization constant. Moreover, if Z pa rticic is in- 
sensitive to the quantum properties of the par- 
ticle, then Xparticic can provide useful informa- 
tion about spectroscopic level ordering Q . Actu- 
ally we shall make use of this result to compare 
the masses before and after the transition and we 
shall not need to be careful about the definition 
of G(t = 0). 

In Fig. 6 the susceptibilities for pions and deltas 
are shown. After the transition the chirality does 
not matter and the masses of parity opposite par- 
ticles coincide. This is another way to see that the 
chiral symmetry is restored. We see again that 
this restoration occurs at a/i ~ 0.4 for (3 — 1.5. 
In Fig. 7 we see the analogous plot displaying the 
same phenomenon for zero density. 
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Figure 4. Polyakov loop II as a function of the 
temperature (written as (3) for zero chemical po- 
tential. 



7. CONCLUDING REMARKS 

We have presented several plots obtained from 
a numerical simulation of QCD with gauge group 
SU(2) at finite temperature and density. Al- 
though our results are still preliminary, we can 
draw a qualitative picture of the phase diagram of 
the theory in the temperature-chemical potential 
plane (we actually substitute the temperature for 
the (3 value, having fixed the temporal size of the 
lattice). There is a region where confinement, chi- 
ral symmetry breaking and a non-zero value for 
the topological susceptibility occurs. This region 
is bounded at /i = by f3 ~ 1.6 and at T = by 
a[i ~ 0.4 ||. We have checked this picture from 
several measurements, all of them giving figures 
in agreement; for instance the position at (3 = 1.5 
where the Polyakov loop raises from zero, the chi- 
ral condensate and the topological susceptibility 
vanish and the particle susceptibilities coalesce is 
always a/i ~ 0.4. 
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Figure 5. Chiral condensate as a function of the 
chemical potential at (3 = 1.5. 
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Figure 6. Particle susceptibility as a function of 
the chemical potential at (3 — 1.5. 
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Figure 7. Particle susceptibility as a function of 
(3 for n = 0. 



